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LONG  TERM  GOALS 

Time  and  space  dependent  radiance  distributions  at  the  sea  surface  are  a  function  of  the  shape  of  the 
incident  distribution  on  the  surface,  modification  by  the  sea  surface  itself  from  topography  and 
transmission  characteristics,  and  alteration  by  the  Inherent  Optical  Properties  (IOPs)  of  the  surface 
ocean.  The  long  term  goal  of  the  proposed  work  is  understanding  this  last  controlling  factor.  With  a 
knowledge  of  the  IOPs,  radiance  fields  can  be  directly  computed  from  the  incident  field  using  the 
equation  of  radiative  transfer,  now  embedded  in  commercially  available  code  (e.g.,  Hydrolight). 

With  the  state  of  current  technology  and  methodologies,  the  primary  obstacles  in  understanding 
subsurface  IOPs  and  their  high-frequency  dynamics  are  a  lack  of  1)  volume  scattering  instrumentation, 
2)  comprehensive  inversion  models  linking  the  IOPs  with  the  causative  bubble,  particulate,  and 
dissolved  matter  in  the  water  (which  in  many  cases  will  require  input  dependent  on  1),  and  3)  suitably 
stable,  non-intrusive  platforms  to  sample  the  subsurface  ocean.  The  first  two  challenges  are  addressed 
in  this  project. 

OBJECTIVE 

There  are  two  overarching  objectives  for  this  project: 

1)  To  develop  an  in-situ  volume  scattering  function  device  measuring  volume  scattering  from  10°  to 
170°  at  10°  intervals  and  sampling  rates  of  1  s'1  or  better  to  sample  the  VSF  in  near-surface  waters;  and 

2)  To  develop  and  refine  IOP  inversion  models  to  resolve  particle  field  characteristics  on  small  spatial 
(cm’s)  and  temporal  (<1  s)  scales  in  near-surface  waters. 
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APPROACH 


Our  design  for  the  VSF  device  is  illustrated  in  Figure  1.  The  device  is  called  MASCOT  (Multi-Angle 
Scattering  Optical  Tool).  The  source  beam  is  a  30  mW  650  nm  laser  diode  expanded  with  a  Gallilean 
2X  beam  expander  to  an  approximately  3mm  X  8  mm  elliptical  shape.  A  wedge  depolarizer  is  used  to 
provide  the  unpolarized  light  needed  for  VSF  determinations.  Seventeen  independent  silicon  diode 
detectors  spaced  in  a  semicircle  10  cm  around  the  sample  volume  measure  volume  scattering  from  10° 
to  170°  at  10°  intervals.  The  total  pathlength  for  all  scattering  measurements  is  20  cm.  Independent 
detectors  allow  resolution  of  the  VSF  without  any  moving  parts  and  time-consuming  scanning. 
Additionally,  each  detector  can  be  optimized  for  its  specific  measurement.  Detector  field-of-views 
(FOVs)  range  from  0.8°  to  5°  for  the  different  detectors,  with  the  narrowest  FOVs  associated  with  the 
detectors  measuring  scattering  in  the  forward  direction.  Using  proprietary  electronics,  a  20  Hz 
sampling  rate  for  all  channels  has  been  achieved  while  maintaining  a  worst  case  signahnoise  of  300:1. 
Relatively  fast  sampling  rates  are  important  in  resolving  VSFs  in  the  highly  dynamic  ocean  subsurface. 

MASCOT  calibration  parameters  have  been  derived  in  solutions  of  microspherical  beads  and  the 
particle  standard  Arizona  Road  Dust  (AZRD;  Powder  Tech.  Inc.).  A  reference  detector  may  be  used  to 
account  for  any  fluctuations  in  source  intensity.  Coincident  measurements  of  beam  attenuation  are 
required  to  correct  for  light  losses  along  the  optical  path  of  the  scattering  measurements  (all 
pathlengths  nominally  -20  cm).  Algorithms  have  also  been  contemplated  that  would  account  for 
attenuation  without  a  separate  direct  measurement.  The  design  in  Fig.  1  is  made  as  small  as  possible 
while  accommodating  all  necessary  hardware,  a  criterion  established  in  part  to  minimize  this  effect. 
Two  MASCOT  prototypes  are  being  developed  for  redundancy  in  case  a  problem  with  a  component 
arises.  Two  systems  allows  for  concurrent  testing  of  different  sensor  components. 

For  the  inversion  modeling,  we  will  extend  the  capabilities  of  existing  models  (Twardowski  et  al. 

2001;  Twardowski  and  Zaneveld,  2004)  by  incorporating  input  from  new  VSF  measurements  and  by 
adding  bubble  particle  populations  (clean  and  coated)  in  the  models.  At  a  minimum,  we  expect  the 
inversion  model  will  give  us  component  concentrations  and  size  distributions,  with  the  components 
being  organic  particles  (living+detrital),  minerals,  and  bubbles.  We  expect  to  co-deploy  the  MASCOT 
and  the  commercially  available  near-forward  VSF  device  LISST  (Sequoia  Inc.)  in  order  to  capture  the 
VSF  with  good  resolution  from  ~0.1  degrees  to  170  degrees. 

RaDyO  field  deployments  are  currently  planned  in  collaboration  with  other  RaDyO  investigators  off 
Scripps  Pier  in  January  2008,  in  Santa  Barbara  Channel  in  September  2008,  and  off  the  Hawaii  islands 
in  May  2009.  For  these  efforts,  we  plan  to  concurrently  deploy  CTD  +  AC9  (+  ACS)  +  MASCOT  + 
LISST  +  ECO  sensors  in  a  real-time  vertical  profiling  system  (Fig.  2).  We  also  expect  to  integrate  two 
additional  sensors  for  these  deployments:  1)  a  high  sampling  rate  fish-eye  lens  radiometer  currently 
being  developed  by  Marlon  Lewis  and  Scott  McLean  of  Satlantic,  and  2)  a  bubble  acoustic  resonator 
provided  by  David  Farmer  and  Svein  Vagle.  For  the  Scripps  Pier  experiment,  a  mobile  deployment 
platform  will  be  configured  so  that  measurements  may  be  made  along  the  length  of  the  pier,  with  the 
expectation  that  higher  concentrations  of  bubbles  will  be  observed  nearing  the  surf  zone.  For 
subsequent  field  work  in  SBC  and  Hawaii,  we  will  deploy  our  sensor  package  off  the  R/V  Kilo  Moana. 
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Figure  1.  Oblique  view  illustration  of  the  MASCOT.  The  VSF  is  resolved  from  10  to 
1 70  degrees  in  10  degree  intervals.  Detectors  are  wedge  shaped  and  arranged  in  a  semi¬ 
circle  on  an  aluminum  frame  to  minimize  reflections  and  perturbation  of  the  water 
sample  in  the  remote  volume  (center  of  semi-circle).  The  source  assembly  includes  a 
30  mW  650  nm  laser  diode,  reference  detector,  beam  expander,  and  wedge  depolarizer. 
Wiring  from  all  the  detector  modules  and  the  source  module  feeds 
to  a  data  handling  unit. 


WORK  COMPLETED 

•  The  first  MASCOT  prototype  was  completed  in  October  2006; 

•  A  summary  of  calibration,  field  deployment,  and  data  processing  activity  is  provided  in  Table  1 
below; 

Table  1.  Timeline  of  MASCOT  activity  in  FY07.  Time  after  May  2007 
has  been  devoted  to  calibration  analyses  and  data  processing. 


OCTOBER 

week  of  the  9th:  delivery  to  WLE 

18-22:  Deployment  in  Pen  River,  ME  (observed  bubbles) 

DECEMBER 

14-15:  Calibration  exercise;  MASCOT  and  an  AC9  in  series  of  suspensions  of 
Arizona  Road  Dust  (AZRD)  and  microspherical  beads  of  1 .992+/-0.025 
um  size. 

JANUARY 

rework  at  WL  HO:  dark  count  offsets  between  air  -  water  improved 

FEBRUARY 

6-7:  Calibration  exercise;  MASCOT  and  an  AC9  in  series  of  suspensions  of 

Arizona  Road  Dust  (AZRD)  and  microspherical  beads  of  1 .992+/-0.025 
um  size;  this  time  with  all  WLE  ECO  BB  sensors  as  well. 

MARCH 

First  2  weeks:  Deployment  in  Hawaii  (SORTIE) 

21-28:  Deployment  in  northwest  Atlantic  (OCVAL:  Endeavor) 

MAY 

4-10:  Deployment  in  Hudson  plume  region  (OCVAL) 
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RESULTS 


Key  results  in  FY07  were  1)  the  MASCOT  prototype  was  fabricated,  tuned  and  tested;  2)  a  viable 
calibration  algorithm  for  the  MASCOT  was  developed;  and  3)  a  substantial  body  of  field  data  was 
collected  from  several  diverse  water  types  around  the  world.  Fig.  2  shows  the  MASCOT  device  in  a 
custom  cage  with  additional  optical  sensors  ready  for  deployment.  Below  is  a  description  of  the 
calibration  method  that  was  developed  and  some  preliminary  field  results. 

To  calibrate  each  of  the  17  raw  signals  to  a  volume  scattering  coefficient,  [5(8)  (nf 1  sr'1),  a  dark  offset 
(DO)  must  be  subtracted  from  the  raw  digital  counts  and  the  result  multiplied  by  a  scaling  factor  (SF). 
Dark  offset  variability  improved  after  the  January  2007  rework.  The  change  in  dark  offset  counts 
collected  in  air  versus  water  ranged  between  10%  and  310%  of  the  approximate  count  levels  measured 
in  the  lab  in  a  c=0.1  nf 1  solution  of  the  particle  standard  Arizona  Road  Dust  (AZRD).  A  good  target  is 
a  few  percent  for  all  channels.  After  the  rework,  the  change  in  dark  offset  counts  between  air  and 
freshwater  ranged  between  6%  and  160%  of  these  count  levels.  For  all  field  deployments,  dark  offsets 
were  measured  in  situ  approximately  daily.  Standard  deviations  of  the  change  in  dark  offsets  in 
seawater  within  any  cruise  period  were  significantly  better  than  the  changes  observed  between  air  and 
water.  As  a  percent  relative  to  approximate  count  levels  measured  in  the  lab  for  c=0. 1  m"1  AZRD, 
standard  deviations  were  33%  and  35%  worst  case  for  the  Hawaii  cruises.  This  improved  substantially 
for  the  North  Atlantic  cruise  (9%  worst  case)  and  the  Hudson  cruise  (1.7%  worst  case).  The  Hudson 
cruise  was  the  only  period  where  isolated  battery  power  was  used  exclusively,  and  this  presumably 
contributed  to  the  excellent  replicability  in  dark  offsets.  A  conclusion  is  that  working  off  battery 
power  with  a  consistent  sensor  suite  configuration  can  provide  sufficiently  reproducible  dark  offsets. 

The  relationship  between  MASCOT  counts  in  any  channel  and  the  particulate  scattering  at  650  nm, 
bp650,  is  expected  to  follow: 

[COUNTS  -  DO]  =  M  *  bp650  *  exp  (-  L  *  (cp650  *  F  +  W)),  (1) 

where  DO  is  the  dark  offset,  M  is  a  linear  slope  that  is  related  to  relative  gain,  L  is  the  pathlength  (0.2 
m),  cp650  is  particulate  attenuation,  F  is  the  fraction  of  particulate  attenuation  that  is  applicable  (i.e., 
we  do  not  want  to  include  the  scattered  light  that  is  collected  within  the  sample  volume  or  detector 
aperture),  and  W  is  the  pure  water  absorption  (water  scattering  is  assumed  negligible)  at  650  nm  (0.34 
m'1).  Fig.  3  shows  an  example  relationship  between  counts  and  cp650  (assumed  approximately  equal 
to  bp650)  for  AZRD  in  one  of  the  lab  experiments.  As  attenuation  increases,  the  rollover  due  to 
attenuation  along  the  path  is  evident.  Detector  saturation  at  -64000  digital  counts  (14  bit)  is  observed 
in  some  channels. 

In  Eq.  1,  there  are  2  unknowns,  M  and  F.  In  theory,  M  and  F  could  be  solved  for  each  channel  using  a 
least-squares  minimization  of  residuals.  There  are  relatively  few  experimental  points,  though,  so  errors 
can  be  expected  to  be  high.  However,  we  know  that  F  should  be  approximately  constant  for  all  the 
channels  (with  the  exception  that  10,  20,  30,  and  40  degree  channels  have  narrower  detector  apertures). 
Thus,  the  normalized  shape  of  each  one  of  the  channel  regressions  should  be  the  same.  We  therefore 
solved  for  F  by  normalizing  all  of  the  channel  regressions  to  area,  essentially  removing  the  slope  (M) 
or  gain  effect.  Using  a  least-squares  minimization  of  the  model  in  Eq.  1  to  all  the  area-normalized  data 
(included  data  from  every  channel)  resulted  in  a  value  of  0.56  for  F.  Thus,  approximately  44%  of  the 
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attenuation  effect  is  being  removed  by  the  collection  of  secondarily  scattered  light  along  the  path  by 
the  sample  volume  and  detector  aperture. 

If  we  invert  Eq.  1  to  solve  for  bp650,  we  obtain: 

bp650  =  [COUNTS  -  DO]  /  M  *  exp  (  L  *  (cp650  *  F  +  W)),  (2) 

If  the  phase  function  of  a  particle  solution  is  known,  P(9),  we  can  then  obtain  (3(9): 

(3(9  =  [COUNTS  -  DO]  *  P(9)  /  M  *  exp  (  L  *  (cp659  *  F  +  W)).  (3) 

The  term  [  P(9)  /  M  ]  is  the  scaling  factor  (SF).  For  microspherical  bead  solutions,  values  of  P(9)  can 
be  detennined  theoretically  using  Mie  theory  and  values  of  M  can  be  determined  experimentally.  With 
F  in  hand,  we  solved  for  M  individually  for  each  channel  using  a  least-squares  minimization  of  the 
model  in  Eq.  1 .  This  was  carried  out  with  a  series  of  microspherical  bead  solutions  for  which  the 
phase  function  may  be  computed  from  theory. 

The  weighting  functions  of  the  17  MASCOT  detectors  were  computed  analytically  using  code  written 
by  Ron  Zaneveld  and  are  shown  in  Fig.  4.  These  weighting  functions  convolved  with  the  P  calculated 
for  the  beads  are  shown  in  Fig.  5.  Scaling  factors  for  each  channel  were  then  computed  by  dividing 
these  P(9)  values  by  the  experimental  M  values. 

Applying  the  calibration  derived  SF  and  F  parameters,  along  with  measurement-specific  DO’s  and 
cp659,  (3(9)  computed  using  Eq.  3  for  the  different  suspensions  from  the  AZRD  experiment  are  shown 
in  Fig.  6.  ECO  VSF  measurements  calibrated  during  the  same  bead  experiment  show  very  good 
agreement.  Applying  the  calibration  derived  SF  and  F  parameters,  along  with  measurement-specific 
DO’s  and  cp659,  calibrated  P(9)  obtained  with  the  MASCOT  in  Hawaii  onboard  the  R/V  KW  are 
plotted  in  Fig.  7.  Concurrent  ECO  VSF  measurements  are  overlaid  with  circles.  The  ECO  VSF 
measurements  show  excellent  agreement  with  the  MASCOT  data.  This  is  clear  when  looking  at  a 
vertical  profile  from  Hawaii  (Fig.  8).  Preliminary  analyses  with  Fournier-Forand  phase  function  fits 
have  also  shown  that  the  measured  VSFs  may  be  represented  well  by  a  particle  population  with  a 
Jungian  type  particle  size  distribution  with  low  bulk  refractive  index,  although  some  overestimation  at 
mid-angles  is  evident.  It  appears  that  the  combination  of  two  or  more  such  populations  may  be 
required  to  achieve  a  precise  fit,  not  an  unexpected  finding. 

IMPACT/APPLICATIONS 

Progress  and  results  represent  important  steps  toward  the  development  of  a  multi-angle,  in-water  VSF 
device.  Knowledge  of  the  Inherent  Optical  Properties  including  the  VSF  can  be  used  to  predict  and 
optimize  the  performance  of  a  host  of  Naval  operations  that  rely  on  divers,  cameras,  laser  imaging 
systems,  and  active  and  passive  remote  sensing  systems.  These  include  mine  countermeasures,  harbor 
security  operations,  debris  field  mapping,  anti-submarine  warfare,  and  search  and  salvage  operations. 
These  measurements  may  also  be  used  in  environmental  monitoring  and  research  applications  for 
determining  particle  composition. 
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Figure  2.  MASCOT  VSF  device  mounted  with  other  optical  sensors  in  a  custom  cage 
designed  to  sample  the  subsurface  domain. 
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Figure  3.  Digital  response  of  MASCOT  scattering  channels  as  a  function 
of  c(650)  measured  with  an  ac-9  in  numerous  suspensions  ofAZRD. 
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Figure  4.  Analytically  computed  weighting  functions  for 
the  1 7  MASCOT  scattering  measurements. 


Figure  5.  Phase  function,  P(  0),  of 1. 992±0. 025  um  beads.  Circles  are 
values  obtained  from  convolving  the  phase  function  with  the  weighting 
functions  for  the  MASCOT  scattering  measurements. 
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Figure  6.  Fully  calibrated  MASCOT  VSFs  obtained  in  suspensions  of  AZRD.  Concurrently 
collected  ECOVSF  VSF  data  are  plotted  using  (+)  symbols. 


Figure  7.  Calibrated  MASCOT  VSFs  obtained  in  Hawaii. 
Concurrently  collected  ECOVSF  data  are  plotted  as  circles. 
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Figure  8.  (A)  A  profile  from  Hawaii  of  the  17  (3(0)  measured  with  the 
MASCOT;  (B)  selected  (3(100),  (3(120),  and  (3(150)  measured  with  the  MASCOT 
and  overlaid  (3(100),  (3(125),  and  P(150)  concurrently  measured  with  an  ECOVSF. 

The  M  ’  and  ‘E’  labels  in  the  legend  refer  to  MASCOT  and  ECOVSF,  respectively. 
Measurements  at  120  and  150  degrees  are  offset  by  0.0002  and  0.0004  m'1  sr,  respectively, 
for  clarity.  Agreement  between  the  MASCOT  and  ECOVSF  measurements  is  excellent. 
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TRANSITIONS 


We  expect  that  our  efforts  in  developing  an  in-water  VSF  device  and  associated  inversion  techniques 
to  better  understand  particle  dynamics  in  natural  waters  will  lead  to  transition  as  operational  tools  for 
the  fleet  and  the  oceanographic  research  community  in  the  future. 

RELATED  PROJECTS 

This  effort  is  related  to  several  ongoing  efforts  to  develop  optical  sensors  and  associated 
biogeochemical  inversion  techniques.  Related  projects  include: 

•  resolving  the  optics  and  dynamics  of  subsurface  bubble  populations  in  the  Southern  Ocean, 

•  developing  novel  harbor  security  monitoring  capabilities  with  Chuck  Trees  and  Jim  Mueller, 

•  developing  improved  vicarious  calibration  and  validation  methods  for  ocean  color  satellite  remote 
sensing, 

•  investigating  the  sources  of  backscattering  in  natural  waters, 

•  developing  tools  for  ocean  observing  systems,  and 

•  developing  a  surfzone  optical  mine  countermeasure  drifter. 
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